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Abstract 



We investigate the collective magnetoplasmon excitations in a quantum 
dot containing finite number of electrons in the high magnetic field limit. We 
consider the electrons in the lowest Landau level and neglect mixing between 
the higher Landau levels. The dispersion relation of these edge modes are 
estimated following the energy weighted sum-rule approach. In this finite 
size system the edge magnetoplasmon modes have different multipolarities 
(or angular momemtum I). Their dependence on the magnetic field and on 
the system size is investigated. With increasing magnetic field, energy of these 
collective modes decreases and in the bulk limit they become gapless. We also 
consider the breathing mode of a dot in the presence of a strong magnetic field, 
and the energy of this mode approaches the cyclotron frequency hto c . 

PACS numbers: 73.20.Dx, 72.15.Rn 
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Recent developments in nanofabrication technology have made it possible to manufac- 
ture electronic systems with reduced dimensionality. Quantum dot is one such nanosystem, 
which is very interesting both experimentally and theoretically. A quantum dot is a two- 
dimensional electron system where the electrons are confined within a finite area by applying 
a gate voltage ]IJ. It is an example of a finite size quantum system where the number of 
electrons can be varied from a few electrons to a few thousand electrons. This nanosystem 
shows many interesting properties in the presence of a magnetic field. One of the interest- 
ing properties in the presence of strong magnetic field is the chiral edge magnetoplasmon 
excitations. The edge modes in the quantum Hall bar have been studied theoretically us- 
ing one dimensional field theory and numerical diagonalisation of the few particle system 
|2|||. The edge modes of a two-dimensional electron system in the presence of positive ion 
background have been studied by several authors ||||. Giovanazzi et al [Q have calculated 
the dispersion relation of the edge modes of the electron gas confined in a jellium disk in 
the long wave-length limit. Recently collective excitations in the finite size nanostructures 
like quantum dot have become very interesting, because many device parameters can be 
controlled by adjusting the gate voltages. The number of electrons in the system can also 
be varied. In particular, the dipole excitation in a quantum dot with parabolic confinement, 
and in the presence of a strong magnetic field has been discussed in ref. pfl. Evidence has 
been found for a strong collective dipole mode in the presence of a magnetic field in far 
infrared spectroscopy experiment |7]||. Recently the collective charge and spin excitation in 
a dot containing about 200 electrons has been observed experimentally ||. Non parabolic 
confinement has been observed in this experimental device, but the exact nature of the con- 
finement is not well understood. The edge magneto-plasmon modes of a quantum dot with 
finite number of electrons have been studied by using magneto-hydrodynamics by several 
authors [0-0. 

The collective excitations in the finite fermion systems like nuclei and metal clusters 
have been extensively studied within random phase approximation(RPA) and using sum- 
rule approach in the last years by several authors [13,15] . This method has also been applied 
to analyze the multipole excitations in the two dimensional quantum dot and antidot systems 
P JT6|Ji7[] . Most of the theoretical works have been done numerically or by using classical 
magneto-hydrodynamics. The main aim of this work is to estimate the collective modes of a 
dot in the presence of a strong magnetic field microscopically, and also to obtain analytical 
results for the dispersion relations of these modes. 

In this paper we study the low-lying multipole excitations and the breathing mode of a 
quantum dot in the presence of a strong magnetic field by using RPA sum-rule approach. We 
consider the electrons in the lowest Landau level, and derive a simple semiclassical energy 
functional for the electrons. By minimising the energy functional we obtain the ground 
state density profile of the electrons. Then we generalise the sum-rules in the presence of a 
magnetic field. Different moments of the RPA strength distribution function are calculated 
by using the ground state density profile. Finally we obtain the analytic expressions for the 
dispersion relation of the low-lying magnetoplasmon modes with different multi-polarities in 
a quantum dot. With increasing magnetic field, the energy of these modes decreases and in 
the bulk limit these modes become gap-less. Neglecting the 1/N corrections, the dispersion 
relation takes a very general form, where only one parameter contains the information about 
the two-body interaction and the shape of the density profile. Also, the 1/N corrections 
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of these modes are calculated. In the strong magnetic field limit the dispersion relation of 
these modes agrees with the dispersion relation obtained from hydrodynamics 1 370f| . We also 



calculate the energy of breathing mode excitation, which is much larger than the energy 
of low-lying edge modes. We obtained the most general form of the dispersion relation of 
the edge excitations in strong magnetic field limit. We also calculated the finite size N 
dependent corrections of the energy of these modes. Analyzing the dispersion relation in a 
strong magnetic field, we have shown the stability of the self-consistent semicircular density 
profile, which shows that this semi circular density is a better ansatz for ground state density 
than commonly used flat density profile. 

The hamiltonian of the two-dimensional electron system in the presence of a magnetic 
field in the perpendicular direction is given by, 



1 . „ . 1 
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+ Enkl-r,|), (1) 

i<j 

where m* is the effective mass, ujq is the frequency of the external parabolic confinement, and 
uj c is the cyclotron frequency Neglecting the two-body interaction, the single particle 
Hamiltonian can be solved exactly. The single-particle energy levels are, 

/ > ) 

E n<m = (2n + \m\ + 1)Q - —u c , (2) 



where = ytOQ + and m is the angular momentum quantum number, and n denotes 
Landau level index. In the presence of a strong magnetic field all electrons occupy the lowest 
Landau level(LLL), energy levels become almost degenarate and they form a band. In the 
lowest Landau level the single particle wave functions are, 



iMO = -^-A=?(r ) m e' r2/2l K (3) 



Vm! k 



o 

where, Iq = \/^q, and z = x + iy is the complex coordinate. In the high magnetic field 
limit the effect of higher Landau levels can be neglected and the many body wave function 
can be constructed out of these lowest Landau level basis. For electrons at the filling factor 
v = 1 the Laughlin wave function can be obtained by constructing the Van-der Monde 
determinant of the above single particle states, 

tp(r u f 2 , ....r n ) = ILi<j{zi - Zj )e-^i Zi * i/21 *. (4) 

In general we can write the ground state many body wave function constructed out of the 
LLL as, 

ip(ri,r 2 ,..ri..) = f(zx, z 2 , ..)e~^ z ^l 2l l. (5) 

We can calculate the energy functional in the lowest Landau level, by calculating the expec- 
tation value of the Hamiltonian with respect to the ground state wave function given above. 
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After doing some algebra and using some properties of the Bargman space [Tj| , we arrive at 
the following energy functional, 

E[ip) = -tku c J \ij\ 2 d 2 r + m*Q(n - u c /2) J r 2 \^d 2 r 

+ \ [ d 2 n [ d 2 r 2 V{\r l -r 2 \)\p{r l ,r 2 )\ 2 d 2 r 2 , (6) 



2 

where, 

\p(R 1 ,R 2 )\ 2 =< n,r 2 , ...ri, ..\Y,^(Ri - ri)S(R 2 - rj)\ri,r 2 , ...fi... > . (7) 

Within Hartree approximation we can write down a simple local density functional for lowest 
Landau level, 

E[p] = J d 2 r[hu c + m*Q(Q - y)r 2 ]p(r) 

+~ / d 2 n f (PripirjVQn-rzDpfo). (8) 



In the case of a flat droplet, the Hartree energy typically goes as ~ N 3 ^ 2 , and the exchange 
correlation energy goes as ~ N, where N is the number of particles in the droplet Jl9 



Hence in the large N limit, the exchange term may be neglected and we obtain the above 
simple looking energy functional. We minimize the free energy E — p,N with respect to p, 
keeping the total number of particles N fixed, and obtain the following integral equation for 
the density, 

m*Q(n - uo c /2)r 2 + j d 2 r'V{\f- r'\)p(r') = p - -hw c . (9) 
For coulomb interaction, V{r) = e 2 /r the solution of the above integral equation is [|1(],|(J, 



P(r)=n 0] /l-(^y, (10) 

where, no = -J^' anc ^ = 8m* 3 n(n-^) • From this simple energy functional we obtain the 
nonperturbative density of the electrons. 

We now consider the edge magnetoplasmon excitations in the dot. Magnetoplasmon 
excitation are the multipole excitations of a quantum dot in the presence of a magnetic 
field. We denote the multipolarity of the operator by the number k. The excitation with 
multipolarity k is generated by the excitation operator F = z k , where z is the complex 
coordinate x + iy. Excitation energies of the above modes can be estimated by sum-rule 
method. Many useful quantities may be calculated from the strength function, that is defined 
as 

S±(E) = E I < n \ F ±\U > ?K E - E n), (11) 

n 

where, E n , and \n > are the excitation energy and excited state respectively, and F + = F, 
and F_ = F^. Various energy- weighted sum rules can be derived through the moments of 
the strength function, and are written below, 
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m$ = \j ' E\S + {E)±S_{E))dE 

= \(< ®\F\H - E ) k F\0 >±< 0\F(H - Eo) k F*\0 >). (12) 

Some useful moments can be written in terms of the commutators of the excitation operator 
F with the many body Hamiltonian H, and they are given below, 

(13) 
(14) 
(15) 
(16) 

where J = [H,F]. In the presence of magnetic field +k and —k collective modes split, and 
the corresponding strength distributions are sharply peaked at these collective frequencies. 
Near the collective excitation energy, we can approximate the strength distribution by delta 
function, S± = a±5(E — E CT ). Since for the multipole modes itlq = 0, we obtain a + = er_. 
Within this approximation the low lying multipole modes can be estimated by using the 
above moments of the strength distribution function and can be written as, 



m = 


i<0|[Ft,F]|0>, 


ml = 


i<0|[Ft,[tf,F]]|0> 


rrq = 


i<0|[Jt,J]|0> 


m t = 


1 -<0\[J\[H,J]]\0>, 



E c zp 



mt 3 m 2 m 2 



where, ± sign correspond to k and —k multipolarities. Similar expression for the excitation 
frequencies of the edge multipole modes can be derived by using the variational principle. 
Given the N electron ground state |0 >, it is possible to find the collective excitation energy 
and the collective state |c >, if one able to find an operator 0\ which satisfies the following 
equation of motion, 

[H,0^]=hu coU Ol (18) 
The state O^|0 > has the excitation energy Tiuj co u. This excitation energy is obtained from, 

_ <0|[Q,[^,Qt]]|0> 

huJcM ~ <o|[o,ot]|o> ■ (19) 

For k — 1 the dipole excitation operator is, 

0^lZ^,~P. + ), (20) 

where p + = p x + ip y . Similarly, we may take the variational ansatz for the higher multipole 
excitations in the following form, 

0^ = F + aJ ; J=[H,F], (21) 
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where a is the variational parameter, F = J2i zf, and J — (~^) Si z f where, D + = 

p+ + %m * 2 u)c z. The collective excitation energy in terms of the the energy weighted sum rules 
is, 



mf + 2am 2 + a 2 mt o 

h^coii = + 2 _ • 22 

2amJ + a z m 2 



After minimising the above expression with respect to the variational parameter a, the 
following expression for the collective excitation frequency is obtained, 



coll 



+ 3 m 2 m 2 

^3- i (^r) 2 + — • (23) 



This agrees with the previous expression in eqn.(17), obtained by the approximate form the 
strength distribution. 

Now we explicitly evaluate the the sum-rules by calculating the commutators. After 
doing some algebra the sum-rules can be written in the following form. The first energy 
weighted moment to* is given by, 

m +(jfe) = < r 2(*-i) >, (24) 

where < ... > denotes the average weighted with the ground state density. The contribution 
for the third moment coming from the kinetic energy of the hamiltonian m^(T) can be 
written as, 

m+(T) = (^)[h 2 k(k - 1) / d 2 rr 2 ( fc - 2 )|Wf 

+2h 2 (k -l)(k-2)< r 2 ^l 2 z > -3h 2 k(k - l)m*cu c < r 2(k ~ 2 % > 
+h 2 m* 2 uj 2 (h 2 + l -k) < r 2 ^ > . (25) 

Similarly the contribution of the external potential term in is, 

m +(V) = ( T ^)km* 2 uj 2 < r 2 ^ > . (26) 

The most important contribution comes from the electron-electron interaction term, and is 
given by, 

m +( ee ) = p^[f d 2 rV 2 V H {r)r 2 ^p{r) + 2{k - 1) / d 2 r^r 2 ^p(r) 

+ I d 2 r I d 2 r'p'(r)e- iW ^-^e iW 'r'^p'(r% (27) 

J J |t" — 7~'| 

where Vff(r) is the Hartree potential. In the presence of a magnetic field there is nonvanishing 
second moment m 2 , 

m - = ^[- h m*uj c k < r 2(k -^ > +2h(k - 1) < r 2(k -Vi z >}, (28) 
m* 2 



where l z is the angular momentum operator with eigenvalue m. We now consider the limit 
in which the magnetic field is large, so that all electrons are in the lowest landau level. 
Neglecting the mixing of higher Landau levels, if we consider the wave function in the 
following form, 

il>(z i ) = f(z 1 ,z 2 ,....z i ,...)e-Z*i«/* a , (29) 

then the above sum-rules can be simplified and written in terms of moments of radius r. In 
the lowest Landau level the sum-rules and can be written as, 

m+(T) + m+(V) = ^[h 2 m* 2 Q 2 (3((k - 1) - ^(|)) 2 + 1) < r 2 ^ > 
-3h 3 m*Q{k - l) 2 {2{k - 2) - fc(^)) < r 2(k ~V > 

+4h\k-l) 2 (k-2) 2 <r 2(k -V >], (30) 

h 2 k 2 

m~ = -\-%m*uj c k < r 2(fe_1) > 

m* 2 

+2h(k - l)(m*n < r 2(fe - 1} > —h(k - 1) < r 2{k - 2) >)]. (31) 
We can write down the expressions for the ratios of moments in the following form, 

^ = 2(t-l)(n-^)- Wc + ^, (32) 
^n^m = + 3( , _ t _ + (33) 

For large number of electrons, if we neglect the terms 0(1/ R 2 ) and the higher order terms, 
then the excitation energies can be written in the most general form, given below, 



Ecoaik) = (k - l)h(Q - p + ^K 2 W + A(e) - ^. (34) 

The coefficients a 1: a 2 , fa depend on the density of the electrons, and the function A(e) 
depends on the nature of the two body interaction and also on the shape of the density profile. 
Now we can estimate the low-lying multipole excitations by using the semi-circular density 
profile. The parameters in the above form of the collective frequency can be evaluated by 
using the semicircular density profile p = n ^Jl — (r/i?) 2 , 



" C)l r(fc + |) 

2 j T(A;)r(3/2) 

a 1 = m(k-l)(2k + l), (36) 



A(e) = 2h 2 n(n - ^)[ \ T 2 / - k], (35) 



a 2 = ~h 2 Q 2 (k - l)(2k + l)(2(k -2)- fc(^)), (37) 



fa = h 2 n 2 (k - l)(k - 2){2k + l)(2k - 1). (38) 

Since the semi-circular density profile is the exact solution of the Hartree equation, the 
dispersion relation of the edge modes obtained by using this density is a non-perturbative 
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result. Further the asymptotic results do not contain the coupling constant e 2 //o, only the 
finite size corrections depend on the coupling constant. This density profile goes to zero at 
the turning point R and the difussive tail is absent in this case. But the difussive tail in the 
exact density distribution can only contribute 1/N corrections to the spectrum of low-lying 
excitations [[RJ. From the energy levels of the non-interacting electrons, we can estimate the 
magnetic field above which all N particles go into the lowest Landau level. This gives the 
condition 2Ktt > (N - - cj c /2)(or co /uj c < l/y/N). If k « N, which is true for the 
low-lying modes, then in the strong field limit we can expand the expression for collective 
modes in terms of A(e)/fo 2 fi 2 , and obtain the following result in strong magnetic field limit 

uj c r(fc)r(3/2) uj c 



This result now be compared with the dispersion relation obtained by, Shikin et al [nj using 
classical hydrodynamics, 



= V^L + ^74-^ c /2, (40) 

r(fc + i/2) 

fci2 »- w o r(ife)r( 3 /2) , l 4i ) 

where Tiuj^ is the same as E c (k) in our notation. If we expand this expression in terms of 
the leading term agrees with the dispersion relation obtained from sum-rule approach. 
For comparison, we may also estimate the edge excitations in a quantum dot by using 
the density profile of the non-interacting electrons at the filling factor v — 1. In the large TV 
limit, the density profile of the system can be approximated by, po9(R — r), where, po = -77 • 
This density profile is very sharp near the edge, and therefore singularities arise in evaluating 
A, but after doin the integrals and then taking the limit (r — > R) at the edge, the singularities 
cancel out and we obtain finite value of the parameter A. The parameter A in this case can 
be derived as, 

^ = 2Kl T^-t^ (42) 

This result agrees with the value obtained by Giovanazzi et al H . In this case the excitation 
energies of various multipole modes vanishes at different values of magnetic field. The 
approximate value of the magnetic field, where excitation energy of the k th mode becomes 
negative is, 

uo )3/2 * ^ 

V/ x/2Mn7kJn±" n 2m- 1 J ' v ; 



where Zo = y m*o> ' -^ ner §y °f the higher multipole modes vanish at lower magnetic field. For 
example the mode with multipolarity k = 10 becomes gapless at a magnetic field ~ 9.65T, 
in a dot with N = 40 and hu>o = bAmeV . Softening of these edge modes indicate the edge 
reconstruction of the dot and the formation of new ground state. This phenomena indicates 
the instability of the flat density. The critical magnetic field where the instability sets in is 
obtained from [|T9"| J^T| , 
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Although A is negative for both the cases, the excitation energy for low-lying modes for 
the semi-circular density profile are positive and asymptotically go to zero, which indicates 
the stability of the ground state. For comparisn, we have shown in Fig.l A for different 
multipolarities, evaluated using two different ground state densities. In Fig. 2, the variation 
of few low-lying modes with magnetic field is shown, using the dispersion law given in 
eqn.(34) and eqn.(35). These modes vanish as ~ 1/B with the increasing magnetic field. 
These low-lying modes are important for the thermodynamics of the system at very low 
temperatures. 

Finally we consider the breathing mode of the dot within the same formalism. The 
breathing mode is excitated by the excitation operater, F = r 2 . The operator J is, 

-ih 

J = [H, F] = ( — -)\p x x + xp x +p y y + yp y ). (45) 



The first moment mf is given by, 



2h 2 

m+ = < r 2 > . (46) 

m* 

The third moment can be written in the following way, 

1 d 2 

™>t = 2-q-i < V\H\V > \v=o, (47) 

where \rj >— e^O >. Now the moment ml^ can be evaluated by using the scaling properties 
of the wave function [ 2"2" |, 

>= e- 2rih2/m *\^(xe- 2rih2/m \ye- 2rih2lm *) > . (48) 
From the above scaling property of the wave function, we obtain, 

mt = 8(^) 2 [< f > + \ m *Q 2 <r 2 >}+ 2(^) 2 E tnt , (49) 

= 8(^) 2 [m*Q 2 <r 2 >]+ 2(^) 2 E mt , (50) 
m m 

where T is the kinetic energy operator and E int is the interaction energy of the system. In 
deriving the second step, we used the properties of the lowest-Landau level wave functions. 
The general expression of the excitation energy of the breathing mode is derived as, 



£7 6 = W4ft a n* + -!L-%-. (51) 
V m* < r 2 > 

Using the semi-circular density for the ground state of the electrons, we obtain the non- 
perturbative result for the breathing mode, 
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E t = 2^1 + 1(1-21). 



(52) 



In the strong field limit this mode approaches the bulk collective mode hoo c . 

To summerise, we have considered the quantum dot in the lowest Landau level. We have 
derived a simple local density functional for the electrons in the LLL within the Hartree ap- 
proximation. For the two-body coulomb interaction a semicircular density profile of the elec- 
trons has been derived. We have estimated the low-lyimg edge multipole modes within the 
sum-rule approach. Since the semi-circular density is exact within Hartree approximation, 
we obtain the non-perturbative results for the dispersion realtions of the collective modes 
by using this density profile. The energy of these collective modes decreases with increasing 
magnetic field. In the strong field limit the most general expression for the dispersion relation 
of the edge modes has been derived and is given by, E(k) = \JVL 2 + A(e) — ^ + (k — 1)uiq^-. 
The parameter A(e) depends on the exact nature of two-body interaction and on the number 
of electrons in the system. We evaluate the parameter A(e) by using exact density, as well 
as by using the density profile of the non interacting electrons at filling factor v — 1. In both 
cases A(e) is negative. In the case of flat density the energy of the low-lying edge modes 
become zero at some magnetic field and the softening of the edge modes indicate the insta- 
bility of the ground state. But for the semi-circular ground state density the low-lying modes 
are positive and asymptotically go to zero, which shows the stability of the self-consistent 
density profile. We have also calculated the breathing mode of a quantum dot, and in the 
strong magnetic field the energy of the breathing mode approaches the cyclotron frequency 
u c . The main results of this paper are, the derivation of the most general dispersion relation 
of the edge modes in a quantum dot in a strong magnetic field. The 1/N corrections of the 
energy of these modes are also obtained. Stability of the self consistent density profile has 
been shown by analysing the dispersion relation of the low-lying edge modes. This result 
shows that the semicircular density is a better ansatz for the ground state density than 
commonly used flat density profile, in the strong magnetic field limit. These edge modes 
are important for the edge excitation and edge reconstruction of quantum dot. They also 
determine the low-temperature thermodynamic properties of the quantum dot. 

I would like to thank M. V. N Murthy for his helpful comments. I would also like to 
thank Tapash Chakraborty and R. Shankar for critical reading of the manuscript. 
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FIGURES 

FIG. 1. The parameter A(e) in units of H 2 ujq, for different multipolar ities k, for a dot with 

2 

N = 50, -pi — = 0.8 and u! c /^>o = 10. The open triangles represent values calculated with 

IQTktJO 

semi-circular density and the solid triangles denote the same parameter for flat density profile. 
FIG. 2. The variation of four lowest multipole modes with magnetic field, in a quantum dot 

with semi-circular density profile. Different values of k denote different multipolarities. 
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